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Abstract. We present line observations of different 
molecular species toward the dense core in L483 and its 
bipolar outflow powered by the Class object IRAS 18148- 
0440. 12 CO(2-l) maps show that the outflow is well col- 
limated and asymmetric, and that its gas is warmer than 
the surrounding cloud by at least a factor of 2. In con- 
trast with the outflows from other Class objects, the 
CHaOH^fc-lfc) lines in L483 do not show strong high ve- 
locity wings, although there is a small (~ 0.3 km s _1 ) 
velocity shift approximately along the outflow direction. 
We do not find evidence for a CH3OH abundance enhance- 
ment in the flow, and the CH3OH lines trace a centrally 
concentrated core which we model, assuming optically 
thin emission, as having a density gradient between r _1 
and r- 1 - 5 for radii between 15" and 100". H 2 CO(2i 2 -l u ) 
lines show strong high- velocity wings with the same distri- 
bution as the outflow, and evidence for a H2CO abundance 
enhancement of a factor of 20 with respect to the ambi- 
ent cloud. At ambient velocities and over the central 40", 
this line presents a strong self absorption and a brighter 
blue peak, a characteristic signature of inward motions. A 
simple analysis of the H2CO line profiles suggests an infall 
rate of 2 x 10~ 6 M Q yr" 1 . 

Combining the results from our observations with pre- 
vious work, we discuss the evolutionary status of IRAS 
18148-0440 and its outflow. The bipolar outflow presents 
some characteristics common to other outflows from Class 
sources, like high degree of collimation, gas heating, and 
H2CO abundance enhancement. However, other charac- 
teristics, like its low velocity, the lack of bright SiO or 
CH3OH outflow wings, and the association with a NIR 
scattering nebula (optically invisible) seem more common 
to outflows from the more evolved Class I sources. As 
IRAS 18148-0440 is a Class object based on its spectral 
energy distribution, we propose that it is more evolved 
than other objects in its class, probably in transition from 
Class to Class I. 



Key words: Stars: formation - ISM: clouds - ISM: in- 
dividual objects: L483 - ISM: jets and outflows - ISM: 
molecules 



1. Introduction 

Class objects are the youngest stellar objects known 



(Andre et al. 1993, 1999). They commonly power bipo- 



lar outflows with extreme characteristics like a very high 
degree of collimation and evidence for shock processing 
of molecular gas even in cases of v ery low stellar lumi- 
nosity (see Bachiller & Tafalla 1999 for a recent review). 
In order for these outflows to evolve into the more quies- 
cent ("standard") outflows associated with Class I sources, 
rapid changes in outflow morphology and kinematics have 
to occur in the few 10 4 yr that Class lasts (Andre et 
al. 1993 ). These changes are most likely associated with 
changes in the source itself, which is undergoing its ma- 
jor phase of assembling via gravitational infall (e.g., Bon- 
temps et al. 1996, Mardones et al. 1997). Understanding 
how these first evolutionary changes of the stellar and out- 
flow life occur is a major challenge to star formation stud- 
ies, and it requires the simultaneous analysis of Class 
objects, their outflows, and their dense gas environments. 
Here we present a molecular line study of the L483 core 
and its outflow powered by IRAS 18148-0440 (IRAS 18148 
hereafter), a system that we find at the end of its Class 
stage, starting its transition to become a Class I object. 

The source IRAS 18148 i n L483 , first identified as an 
embedded object by Parker ( |l988a ), is one of the reddest 
low-mass sources known (Ladd et al. 1991a, 1991b), and is 
located toward the Aquila Rift, at a most likely distance of 



200 pc (Dam e fc T haddeus |1985| ). Ladd et al. (|991aj) and 
Fuller et al. (1995) estimate a source bolometric tempera- 
ture (in the sense of Myers & Ladd |l993| ) of 50-60 K, and 
using the flux compilation by Fuller et al. (1995) (their 
Fi g 4) we estimate a L smm /Lb \ > 0.9 % (also, Fuller et 
al. 1995 fit the spectral energy distribution with a single- 
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temperature dust model at 40 K, and from their 1.1mm 
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flux, we estimate ibol/ii.imm < 2.5 x 10 4 ). These num- 
bers suggest that IRAS 18148 is a Class object (Andre 
et al. 1993, 1999, Chen et al. 1995), as already proposed 
by Fuller et al. ( |1995[ ) and Fuller & Wooten ( |2000[ ), al- 
though it is less extreme than the prototype Class source 
VLA1623 (Andre et al. [l990| Andre et al. p993| ). IRAS 
18148 has a luminosity of about 10 L Q and drives a wcll- 
collimated bi polar CO outflow (Par ker et al. |1988b| , [1991 , 



Fulle r et al. |1995| , Bontemps et al. |1996| , Hatchell et al 



199S ), and is associ ated with a variable H2O maser (Xi- 
Turner 1995 ) a nd sh ocked H2 emission (Fuller et 



ang 

al. 1995 , Buckle et al. 199£ ). NIR imaging of the source 



vicinity shows a well-defined, parabolic reflection nebula, 
which is optically invisible and coincides with the blue 
lobe of the CO outflow (Hodapp |1994| , Fulle r et a l. |l995| ). 
Ammonia observations by Goodman et al. (1993), Fuller 
& M yers (|1993| ), Anglada et al. ( |1997| ), and Fuller & Woot- 
ten ( [2000 ) reveal that the L483 core is centrally concen- 
trated, has a strong velocity gradient across it, and a gas 
kinetic temperature of about 10 K. H2CO and CS spectra 
toward the central source present strong self absorption 
with lines having brighter blue peak, a signature of infall 
motions (Myers et al. 1995| , Mardones et al. 1997). 

The combination in L483 of Class characteristics, like 
a low Tboi and infall asymmetry, together with the pres- 
ence of a bright NIR nebula, indicative of partial core dis- 
ruption, makes this source an interesting object to study 
the early evolution of a very young stellar object. To carry 
out such a study, we have observed L483 in tracers sensi- 
tive to different aspects early stellar life, like the outflow 
(CO, section 3.1), the dense core and possible chemical 
outflow anomalies (CH3OH, section 3.2), and infall and 
shock chemistry (H2CO, 3.3). From the combination of 
these observations, we propose that the central source in 
L483 has already started its transition toward Class I, 
and that the outflow has lost part of the chemical richness 
characteristic of Class flows (section 4). 

2. Observations 

We observed L483 with the IRAM 30m telescope dur- 
ing several sessions in 1994 September, 1995 September, 
November, and 1996 June. Different receiver configura- 
tions were used to map the core in 12 CO(2-l) [230.53799 
GHz], H 2 CO(2i2-ln) [140.839518 GHz], and CH 3 OH(2 fc - 
lfc) [96.741420 GHz], and to observe selected positions 
in 13 CO(1-0), C 18 O(l-0), C 17 O(l-0), Hi. 3 CO(2i2-ln), 
H 2 CO(3i2-2n), and SiO(2-l). Most observations were 
done in position switching mode (PSW) in order to ob- 
tain flat baselines. After searching for a clean off position, 
we settled with (—600", 300) with respect to our map 
center (ai 950 = 18 h 14 m 50!6, S 195 „ = -4°40'49"0, position 
of IRAS 18148-0440). This position lies outside the opti- 
cal obscuration associated with L483 and seems free from 
dense gas tracer emission, although it has a weak 12 CO(2- 
1) line at a level of 2 K between Vlsr = 6 — 9 km s _1 



(outside the L483 range). A frequency switched (FSW) 
12 CO(2-l) spectrum of this position was taken so it could 
be added to the data if needed, and a test was made by 
adding the off spectrum to a PSW spectrum from the ori- 
gin and comparing the result with a FSW spectrum; the 
two were indistinguishable. The data shown in this paper, 
except when indicated, correspond to PSW observations 
without addition of the off position. 

The backend was an autocorrelator split into differ- 
ent windows with resolutions ranging from 0.1 km s _1 for 
12 CO(2-l) to 0.03 km s^ 1 for CH 3 OH(2 fe -l fe ). The tele- 
scope pointing was corrected frequently by observing con- 
tinuum sources and is expected to be accurate within 3". 
The T£ scale of the telescope was converted into T m b using 
the main beam efficiencies recommended by Wild ( 1995 ) . 
The full width at half maximum of the telescope beam 
varies linearly with wavelength, and ranges from 25" for 
CH 3 OH(2 fc -l fe ) to 11" for 12 CO(2-l). Spectra were taken 
with 10" spacing, which is slightly less than one beam at 
the highest frequency. 



3. Results 

3.1. CO data 

3.1.1. Ambient CO emission 

In this section we derive the basic parameters of the gas 
along the line of sight towards the L483 core center us- 
ing a series of CO isotopomer lines. As Fig. [l] shows, the 
12 CO(2-l) line is heavily self absorbed between Vlsr 4.5 
and 6.5 km s , with a T m b of 4-5 K Given that the gas 
between these velocities is optically thick, its excitation 
temperature T cx = has to be 8.5-9.5 K, a value close to 
the kinetic temperature of th e den se gas (10 K, fro m NH 3 
and HC 3 N, Fuller & Myers |1993| , Anglada et al. |1997| ). 
The gas responsible for the 12 CO(2-l) absorption repre- 
sents low density gas in the outer layers of the cloud, far 
from the dense core, so the above result suggests that all 
the ambient gas along the line of sight has an almost con- 
stant temperature of 10 K (note that the outflow gas is 
warmer, see below). 

In contrast with the self absorbed 12 CO(2-l) emis- 
sion, C 18 O(l-0) and C 17 O(l-0) are optically thin. This 
is not evident because of the non-Gaussian shape of the 
spectra (Fig. [I]), but can be proved using the following 
simple model: we create a C 17 O(l-0) spectrum by adding 
three replicas of the C 18 O(l-0) line each shifted in ve- 
locity by the proper hyperfine splitting and weighted by 
the optically thin relative intensity. The result, shown 
in lighter shade below the C 17 O(l-0) spectrum, matches 
very well the observations, indicating that both C 18 0(1- 
0) and C 17 O(l-0) are thin. The non Gaussian shape of 
these lines has therefore to result from velocity structure 
along the line of sight (see section 3.2 for further details). 

With the excitation temperature of 10 K and the fact 
that C 17 O(l-0) is thin, we can derive the core central 
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Fig. 1. Spectra of CO isotopomers toward IRAS 18148- 
0440, at the center of the L483 core (ai 950 = 18 h 14 m 50!6, 
5 W50 = -4°40'49"0). For 12 CO(2-l), the spectrum toward 
the non outflow position (0, —40") is also shown in lighter 
shade to better illustrate the deep self absorption and out- 
flow wings toward the central position. Both 12 CO(2-l) 
spectra have been corrected for emission in the off posi- 
tion. Below the C 17 O(l-0) spectrum, and in lighter shade, 
is shown a model using three replicas of the C 18 O(l-0) 
line with the relative weights expected for optically thin 
emission. Note the good agreement, which indicates a neg- 
ligible optical depth in C 17 O(l-0). 



H2 column density. We integrate the C 17 O(l-0) emission 
in velocity and assume local thermodynamic equilibrium, 
estimating a N(C 17 0) of 1.5 xlO 15 cm~ 2 . For a standard 
C 17 abundance of 4.7x10" 



3 xlO 22 cm- 

To finish this section, we note the presence in Fig. [j] 
of CO emission outside the ambient cloud range (Vlsr, = 
4.5-6.5 km s _1 ). Part of this emission comes from outflow 
gas (discussed in section 3.1.2), but other part must come 
from additional clouds along the line of sight. This is the 
case for two features at Vlsr = 3 and 8 km s _1 , mostly 
seen in 13 CO, as they are spread over an area larger than 6 
arcminutes. It must also be the case for some contribution 
between 2.5 and 4.5 km s _1 , which is very prominent in 
the 12 CO spectra outside the outflow range (see (0,-40) 



(Frerking et al. |1982[ Wilson 
& Rood 199J), this value implies an H2 column density of 



spectrum in Fig. [l]). As we will see below, this low velocity 
emission limits our study of the outflow gas, to which we 
now turn our attention. 

3.1.2. The outflow 

Fig. H presents our CO(2-l) map of the L483 outflow (see 



also Fig. A.2). As previous maps (Parker et al. 11991 
Fuller et al . [1995| , Bontemps et al. fll996| ), and Hatchcll 



et al. 1999| ), it shows that the accelerated CO emission 



extends E-W with IRAS 18148 at the center. The red CO 
lies to the east of IRAS 18148 while the blue CO lies to 
the west, although there is some red gas near the western 
tip of the blue lobe. This anomalous red gas coincides with 



a region of shocked H2 emission (Fuller et al. 1995), and 
probably represents outflow material with an enhanced 
turbulent component because of the shock. 

A comparison of the CO outflow with the K' (2.1 /xm) 
image of Hodapp ( 1994[ ) (Fig. ^) shows that the blue CO 
delineates very closely the reflection nebula and presents a 
relative maximum toward its center. In contrast, the red 
(eastern) CO is weak and highly collimated near IRAS 
18148 and reaches a maximum 0.04 pc away from it. This 
east-west asymmetry of the CO outflow suggests that a 
similar asymmetry occurs in the underlying nebula, which 
would therefore have an intrinsically brighter western side. 

As the 12 CO (2-1) spectrum in Fig. [i] shows, the out- 
flow gas is warmer than the ambient core, or otherwise its 
emission would not be brighter than the self absorption by 
10 K gas. The strongest CO wings imply excitation tem- 
peratures of 20 K (twice the ambient kinetic temperature) , 
and this is a lower limit because the CO emission may not 
be thick and thermalized. In fact, Hatchell et al. ( 1999| ) 
have argued, from a comparison of CO(4-3) and CO(2-l) 
lines, that temperatures up to 50 K could be present in 
the CO outflow. Given that we do not find brightness tem- 
peratures larger than 20 K, we take this value as a lower 
limit for our further calculations of the outflow energetics 
(see below) and molecular abundances (sections 3.2 and 
3.3). 

With the above value for the CO(2-l) T ex , we estimate 
the outflow energetics. This requires special care because 
the outflow is rather slow and the CO spectra are contam- 
inated at low velocities, so in the Appendix we present 
a method to deal with the contaminating emission. Of 
course, we cannot correct for the outflow gas hidden by 
the self absorption, so all CO emission in the central 3 
km s _1 (about Vlsr = 5.5 km s" 1 ) will be ignored, and 
our estimate will represent a lower limit. As the highest 
outflow velocities, we take Vlsr —3 and 14 km s _1 , be- 
cause no faster outflow emission is found in the CO(2-l) 
velocity maps. Assuming a CO abundance of 8.5 x 10~ 5 
(Frerking et al. |1982| ), we derive an outflow mass of 0.01 
M , a momentum of 0.03 M Q km s , and a kinetic energy 
of 2 x 10 42 erg. From a total outflow length of 160" and a 
total velocity extent of 15 km s , we derive a kinematical 
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Fig. 2. High velocity 12 CO(2-l) emission superposed to the K' (2.2 /im) image of Hodapp (1994) illustrating the 
excellent agreement between the outflow distribution and the scattering nebula. Solid contours represent blue emission 
integrated from Vlsr to 2 km s _1 and dashed contours represent red emission integrated from Vlsr 8 to 10 km s" 1 . 
For both lobes, contours are at 3, 6, ... K km s" 1 . Astrometry for the K' image has been done by finding common 
stars in the DSS image, and the estimated positional accuracy is less than 2". The K' image has been convolved with 
a Gaussian beam to enhance the diffuse emission from the red lobe. The star sign marks the position of IRAS 18148 
and the scale bar is for an assumed distance of 200 pc. 
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Fig. 3. CH30H(2fc-lfc) spectrum toward the center of the 
L483 core showing the different k components (see Fig. [I] 
for coordinates). 



time of 10 4 yr. These values are in reasonable agreement 



with those from from Parker et al. fll99l| ). 



3.2. CH 3 OH and SiO data 

Figure [| presents a CH30H(2fc-lfc) spectrum toward the 
core center with labels indicating the different k compo- 
nents. Among the E-type lines, the k=l component is not 
detected, while the the integrated intensity ratio between 



the 2_ 1 -1„ 1 and 2 -lo lines is approximately 7, very sim- 
ilar to the ratio found by Turner ( |l998|) for TMC-1 and 
L183. The relatively high intensity of the lines (about 
2.5 K at the peak) implies an excitation temperature of at 
least 5.5 K, while a n LTE rotation diagram analysis (e.g. 
Menten et al. 



indicates temperatures of 4 K at most. 
This suggests that non LTE conditions may apply, and 
that for the E-type CH3OH, the lower lying k=— 1 ladder 
has a higher T cx than the higher fc=0 ladder. The excita- 
tion conditions are probably constant over the core, as an 
average over all (75) positions outside the central 20" gives 
a spectrum with the same intensity ratios among all the E 
and A components as in the spectrum shown in Fig. we 
will use below this apparently constant conditions to infer 
a lower limit to the density gradient in the core. From the 
integrated intensities in the central spectrum, we estimate 
a CH3OH column density of about 7 x 10 13 cm" 2 , which 
together with our estimated H2 column density implies a 
CH3OH abundance of 2 x 10 -9 , which is very close to the 
abundances estimated for other dark clouds (Friberg et al. 
198$ Bachiller et al. |1995| , Turner |1998| ). 

The spatial distribution of the CH30H(2fe -lj.) inte- 
grated intensity, shown in Fig. ^a, is rather round and 
centrally concentrated toward IRAS 18148. This is also 
the case with other dense gas tracers like HC3N, NH3, 
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Fig. 4. a Map of CH30H(2 fe— lfe) intensity integrated over the 2o-lo A and 2_x~ l-i E components. Levels go from 0.3 
to 2.7 K km s -1 in steps of 0.3 K km s^ 1 . b Radial average of integrated intensity (filled squares represent 2_i— l_i E 
and open squares represent 2o~lo A) with optically thin, LTE models for different density laws convolved with a 25" 
Gaussian (solid lines), c Position- velocity diagram along R.A. for AS = (i.e., along the outflow axis). Contours go 
from 0.5 to 2.5 K in steps of 0.5 K. Note the velocity shift near Aa = 30". Offsets with respect to IRAS 18148 (see 
Fig. @) 



and N 2 H+ ( Fuller & Myers [1993[ Anglada et al. [1997 , 
Caselli et al. |2000 ), and is probably a sign of the extreme 
youth of the embedded star, which has not had time to 
perturb the bulk of the parental core (despite carving the 
nebula). The central concentration of the emission rein- 
forces our interpretation that the emission is mostly opti- 
cally thin, as otherwise we would expect a flat distribution 
(our unpublished NH3 data indicate a constant gas kinetic 
temperature across the core). 

Fuller et al. (1995) have inferred an r -2 density gradi- 



ent within 30" of the IRAS source based on the intensity 
contrast between the two sides of the IR reflection neb- 
ula, but their calculation depends critically on assuming 
that the nebula is symmetric, which we have seen is prob- 
ably not the case. More recently, Fuller & Wootten (2000) 
have proposed that the r -2 profile continues to larger radii 
(100") using a model for the NH3(1, 1) emission, although 
this result could be affected by interferometer missing flux. 
Here we use the CH3OH emission to derive an indepen- 
dent estimate, assuming that this emission is thin and the 
excitation temperature is constant. In Fig. ^b, we show a 
radial average of the emission, together with the results 
from our model for three density power laws: r^ 1 , r~ 15 , 
and r~ 2 (results convolved with a 25" Gaussian beam). 
As the figure shows, the r~ 2 density law is too steep, 
and the observations are better fit between R = 15" and 
100" using a r _1 profile with probably some steepening 
for R > 50". 

In the optically thin limit, the brightness radial pro- 
files of Fig. |] are proportional to column density radial 
profiles, so we can use them to estimate the core mass by 
normalizing them to the central H2 column density and in- 
tegrating them radially. In this way, we derive a core mass 
between 5 M Q (r~ 15 density profile) and 10 M Q (r _1 den- 
sity profile). These values are in goo d agre ement with the 
ammonia result from Anglada et al. 1997 , and also agree 
with a virial estimate using the average CH3OH line width 



over the core (0.64 km s 1 ), which gives 7 and 8 M© for 
r~ and density profiles, respectively. 

Although the CH3OH lines are relatively narrow com- 
pared with other dense gas tracers like H2CO, they are 
systematically asymmetric and change velocity across the 
core. This is illustrated in Fig. ||c with a position-velocity 
diagram along the east-west axis of the core, which shows 
that the line center velocity shifts abruptly by about 0.3 
km s _1 near the position of IRAS 18148 (origin of offset 
coordinates); lines toward Aa < 20" are brighter, bluer, 
and narrower than toward the east. These two velocity 
components coincide with the C 18 (and C 17 0) peaks we 
have found in section 3.1.1, and and can also be seen in 



our N 2 H+(l-0) data. Fuller & Myers Q1993Q , with low res- 
olution observations, have reported that the spectra from 
L483 seem to have two velocity components, while Good- 
man et al. 1993 found a systematic velocity gradient across 
the core, all in the same direction as the velocity change 
we find. The origin of this behavior is not clear, but we 
notice that the sense of the velocity change agrees with the 
sense of the bipolar outflow. It is therefore possible that 
it arises from the acceleration of dense gas by the out- 
flow, since the position velocity diag ram i s similar to that 
of L1228 in C 3 H 2 (Tafalla & Myers |l997| ), where outflow 
acceleration is the cause of a similar velocity shift. 

Independently of the origin of the velocity shift, it is 
clear that the CH3OH spectra in L483 do not show the 
prominent high- velocity wings seen in some low mass out- 



flows like NGC 1333-I RAS 2 (Sandell et al. |1994[ ) and 
L1157 (Bachiller at al. 1995). The wings in these systems 
arise from large abundance enhancements of CH3OH, 
which in the best studied case of L1157, amounts to a 
factor of 400 (Bachiller et al. |1995[ Bachiller & Perez- 



Gutierrez 1997). For L483, the non detection of CH 3 OH 



at the velocities with CO outflow implies that any pos- 
sible CH3OH abundance enhancement cannot be larger 
than around 10, which is at least 40 times lower than in 
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Fig. 5. H 2 CO(2i 2 -lii) spectra toward selected posi- 
tions of the L483 core. For the top panel, the weaker 
H2 3 CO(2i2-ln) is also shown. Offset are in arcsec with 
respect to the position of IRAS 18148 (see Fig. 0). 



L1157. Higher signal-to-noise data can probably lower this 
limit by a significant amount. 

Another tracer with prominent wings toward certain 



outflows is SiO (e.g., Bachillcr & Perez- Gutierrez 1997), 
and three L483 positions were observed in the 2-1 line of 
this molecule (origin, 30" E, and 30" W). Our non detec- 
tions (with limits of the order of 0.1 K km s _1 ) imply SiO 
abundances lower than 10 -11 and 8 x 10~ 10 for the am- 
bient gas and the outflow, respectively (numbers derived 
using an large velocity gradient analysis). Although the 
non detections make impossible to estimate an abundance 
enhancement in the outflow, comparing our outflow limit 
with the SiO abundance found in L1157 (10~ 7 , Mikami 



et al. 1992 , Bachillcr & Perez-Gutierrez 1997 ), we esti- 
mate that any possible enhancement in L483 is at least 
100 times smaller than in L1157. 



3.3. H 2 CO Data 

Our other dense gas tracer, H2CO, presents very different 
lines than CH 3 OH. Fig. || shows a series of H 2 CO(2i 2 -lii) 
spectra for different core positions illustrating the variety 
of profiles. At ambient velocities (Vlsr ~ 5-5 km/s), the 
lines are strongly self absorbed over most of the core, as 
can be checked by comparing the double-peaked main iso- 



tope with the single-peaked rare H^CO line toward the 
central position (Fig. || top panel). This self absorption 
is slightly red shifted with respect to the emitting gas, so 
the resulting spectrum has a brighter blue peak, a fact 



already noticed by Myers et al. (1995) and Mardones et 



al. (1997), who have presented preliminary versions of the 
H2CO spectrum towards the core center. Such red-shifted 
self absorptions are characteristic signatures of inward mo- 
tions (e.g, Leung & Brown 1977), and their presence in 



L483 makes this core one of the best infall candidates 
known. A detailed study of the spatial distribution of the 
self absorption and its interpretation in terms of infa ll mo- 
tions will be presented elsewhere (Mardones et al. 2000, 



see also Mardones 1998 ), so here we limit ourselves to 
briefly comment on this feature. 

Double-peaked H2CO spectra abound in the core and 
are more prominent along the outflow axis, although a 
map of the spectral asymmetry parameter Sv (defined as 
the difference between the thick and thin line peaks nor- 
malized to the thin line width, see Mardones et al. 1997| ) 
shows that the blue asymmetry is stronger perpendicular 
to the flow. Overall, despite the presence of high-velocity 
blue and red wings in the H2CO lines, the ambient self 
absorption is mostly red shifted and an average spectrum 
over the central core is clearly asymmetric in the sense of 
infall. To estimate the global infall rate in L483, we first 
determine the infall radius from the extent of the H2CO 
spectra with brighter blue peak, which we measure from 
the data as 0.02 pc. Then, we use the simple 2-layer model 
of Myers et al. 



(1996) 



which allows to derive an infall ve- 
locity from the contrast between the blue and red peaks of 
a self absorbed line profile knowing the intrinsic line width 
from a thin tracer (see their Eq. 9). To do this, we take the 
average H2CO spectrum inside the infall radius, and use 
as thin tracer the average CH3OH spectrum over the same 
area. In this way, we derive a mean infall speed of 0.02 km 
s _1 , which is clearly subsonic. Finally, we derive a mean 
density inside the infall radius using the power-law density 
models derived in our CH3OH analysis, which give a mean 
density of 3.3 x 10 5 cm -3 for a r _1 profile and 2.8 x 10 5 
cm~ 3 for a r -15 profile. Averaging the above densities to 
3 x 10 5 cm -3 , we derive a mass infall rate of 2 x 10 -6 M© 
yr _1 (sec Mardones et al. 200C for further details). 

At velocities outside the ambient range, the H2CO 
spectra present strong wings that change with position 
following the distribution of accelerated gas in the bipolar 
outflow. The spectra at (30", 0) and (-10", 0) in Fig. | 
illustrate this effect, which implies that the bipolar out- 
flow is accelerating part of the dense core gas, probably 
shocking it and altering its chemical composition (see be- 
low) . To better illustrate the dense-gas acceleration traced 
by H2CO, we present in Figure ^ velocity maps for the 
blue, ambient, and red regimes. The blue and red maps 
agree very well with equivalent outflow maps from CO (see 
Fig. ||), especially for the blue lobe. There, both emissions 
present two maxima, one associated with the reflection 
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Fig. 6. H2CO(2i2-lii) velocity maps showing blue, ambient, and red emission (left, center, and right panels). Note 
how the H2CO distribution follows that of the CO outflow (Fig. |2|). In all panels, the the first level and the level step 
is 0.2 K km s^ 1 . Offsets as in Fig. 4a. 
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Fig. 7. Position- velocity diagram of the H2CO(2i2-ln) 
emission along R.A. for AS = (i.e., along the outflow 
axis). Note the prominent outflow wings increasing in ve- 
locity with position; compare with Fig. 4c. First level and 
step are 0.2 K. 



nebula, and the other with the region of strong H2 emis- 
sion at the end of the lobe. Even more, both emissions 
present an "anomalous" red peak towards the position of 
bright H 2 emission, again suggesting this gas is shock re- 
lated. These similarities suggest that CO and H 2 CO, de- 
spite their different dipole moments, are tracing the same 
(or very closely connected) gas. 

As a further illustration of the importance of the out- 
flow acceleration in the H2CO emission, we show in Fig. [7] 
a position- velocity diagram along the the Right Ascension 
axis (i.e., approximately parallel to the outflow). At am- 
bient velocities, the strong self absorption can be easily 
seen in the form of two bright emission peaks at about 5.2 
and 6 km s _1 . At high velocities, the outflow wings form 
two triangular-shaped extensions which illustrate how the 
emission terminal velocity increases almost linearly with 
distance from the outflow source (at Aa = 0); such lin- 
ear increases are common in the CO emission from out- 



flows (e.g., Meyers-Rice & Lada 1991). We can also see 
in the position-velocity diagram that towards the west 
(Aa ~ —50"), both the wing and the ambient emission 
drop simultaneously. This suggests that the high velocity 
gas is breaking through the dense core, something also sug- 
gested by the presence of bright H2 emission. Towards the 
east, the ambient gas is more extended than the outflow 
wing, so it seems the outflow is bounded in this direction. 

We finish studying the H2CO abundance in the outflow 
and comparing it to the abundance in the ambient cloud. 
As the H2CO(2i2-ln) line is self absorbed at ambient ve- 
locities, we use in this range the H3, 3 CO isotopomer, which 
is very likely thin (see Fig. ||). In this way, using an LTE 
analysis with T c x = 5 - 10 K and a 12 C/ 13 C ratio of 77 
(Wilson & Rood 1994), we derive a H2CO ambient column 
density of 5 x 10 la cm -2 for the central position (a simi- 
lar analysis using H2CO(3i2-2n) gives the same number). 
This value implies an ambient abundance of 1.5 x 10~ 9 . 
For the outflow regime, we use the H2CO(2i2-ln) line and 
compare it with the outflow CO emission. Applying the 
same LTE analysis as before, we estimate outflow abun- 
dances at (-10", 0) (blue) and (30", 0) (red) of about 
3 x 10~ 8 , suggesting an H2CO abundance enhancement of 
a factor of 20. This number is slightly smaller that the fac- 
tor of 60-80 found in L1157 (Bachiller & Perez- Gutierrez 



1997), but is significantly larger than our limit for the 
CH3OH enhancement. Thus, in contrast with L1157 for 
which CH3OH is enhanced by almost an order of magni- 
tude more than H2CO, L483 is richer in H2CO. 

4. Evolutionary status of IRAS 18148-0440 

As mentioned in the Introduction, IRAS 18148 belongs 



to Class (Andre et al. [1993] Andre et al. |l999j) due to 
its red spectral energy distribution, although it is less ex- 
treme than objects like VLA1623. IRAS 18148 has ad- 
ditional Class characteristics like the presence of infall 



asymmetry in line spectra (Mardones et al. 1997). After 



analyzing its molecular environment in previous sections, 
we now compare this object and its outflow with other 
Class sources, and to do this, we present in Table 1 a 
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Table 1. Class sources and their outflows 









Outflow 


Chem. Enhancement 




Name 


Lboi 


Tboi 


t /-b nc 
•max *-* 


SiO 


CH 3 OH 


H 2 CO 


Ref 




L 


K 












L1448-C 


9 


56 


70 Y 


> 10 5 


7 


? 


1,2,3,4 


L1157 


11 


62 


12 N 


> 10 5 


400 


75 


2,5,6,7 


BHR 71 


9 


56 


10 N 


200 


30 


? 


2,8,9 


L483 


10 


50 


7.5 N 


ND d 


< 10 


20 


2, 10 



With respect to ambient cloud. 

Outflow terminal velocity from CO in km s . 

Presence of CO "bullets." 

Not detected. Enhancement at least 100 times smaller than 
in L1157. 



References: (1) Bachille r et al . 1991a, (2) Mardones et al 
(3) Bachiller e t al. | l990| 



Umemoto et al. 



1992 



(7) Avery & Chiao 



1997 



16 



199' 



et al. [l998| , (10) Fuller et al. |l995 



(4) Bachiller et al. |l9911j, (5j 
Bachiller & Perez -Guti errez [1997 , 
(8) Bour ke et al. |l997|, (9) Garay 



summary of properties of different Class sources with 
similar luminosity (~ 10L Q ) for which enough molecular 
data are available. 

As Table 1 shows, L483 is very different kincmatically 
from L1448-C, a source characterized by an extremely col- 
limated and fast bipolar CO outflow with disc rete h igh 
velocity components ( "bu llets" ) (Bachiller et al. 1990| , see 



also Barsony et al. 1998 



1999, Eisloffel 



O 'Linger et al. 
200C, and the similar object VLA1623 found by Andre et 
al. 1990). Despite our limited knowledge of the chemical 
composition of the L1448-C outflow, Table 1 shows that 
L483 also differs from L1448-C in its chemical enhance- 
ments, again in the sense of L1448-C being more extreme. 
More similarities are found between L483 and LI 157 



an d BH R 71 (L1157: U memo to et al. |1992| , Mikamiet 
al. [1992|, Bachiller et al. |1995[ Tafalla & Bachiller |1995 , 

Bachiller & Perez- 



Zhang et al. [1995| , Avery & C hiao |1996 



Gutierrez 1997 Gueth et al. 1998", Umemoto et al. 1999 



would represent a rather late stage of a Class source, 
and the objects in Table 1 would be ordered by increasing 
age from top to bottom. Unfortunately, no independent 
stellar clock exists yet to order by age the different Class 
objects, and as Table 1 shows, the bolometric temper- 
ature Tboi cannot distinguish between extremely young 
sources, probably due to the lack of mid and far IR pho- 
tometry. There is a reason, however, to suspect that IRAS 
18148 is more evolved than other Class sources, and that 
the above evolution s cenar io is correct: L483 has a bright 
NIR nebula (Hodapp |l994| , Fuller et al. |1995| ) in co ntrast 
with L1448-C an d L11 57 (Bally et al. |l993| , Hodapp |1994 , 



Davis & Eisloffel |1995) (bu t note BHR 71 also has a NIR 



nebula, Bourke et al. 1997 ). The presence of such a neb- 
ula suggests that the L483 outflow has been accelerating 
ambient material for longer than the outflows from the 
other objects in Table 1, and that IRAS 18148 is more 
advanced in its transition to become a visible object (e.g., 
Shu et al. 1987). Other C lass object with a NIR neb- 
ula is L1527 (Eiroa et al. 1994), so if the above scenario 



is correct, one would expect to find in this object a weak 
abundance enhancement of the molecules shown in Table 
1. Observations of this object should be done to test this 
point. 

The status of L483 as a somewhat evolved Class 
object is also consistent with its relatively low value 
L smm /L ho i ( > 0.9 %, compare with the 10 % of VLA1623, 
Andre et al. 1993) and with the work of Bontemps et al. 



Zh ang e t al. ^000[ BHR 71: Bourke et al. |1997| , Garay et Andre et al 



al. 1998 ). These outflows have comparable kinematics in 
the sense of having lower velocities and lacking "bullets" 
(also their collimation is similar), but differ by the amount 
of the chemical enhancement in the sense that the num- 
bers for L483 are systematically smaller. Although the un- 
certainties in the chemical estimates are rather large, the 
differences shown in Table 1 are too extreme to be due to 
observational error. This means that among "chemically 
active" outflows, L483 is a weak case. 

Bachiller & Tafalla 199E have argued that outflows 



1996. Applying the factor of 10 correction these authors 
apply to their CO data (see their Eq. 2), we estimate for 
L483 a momentum flux d la Bontemps et al. of 3 x 10 -5 M© 
km s _1 yr . This value is almost a factor of 2 lower than 
the mean momentum flux of Class objects, but still 8 
times larger than the average number for Class I sources. 
L483, again, appears as a Class object already evolving 
toward Class I. 

If the Class encompasses sources as diverse as those 
powering the L1448-C and L483 outflows, outflow evolu- 
tion should occur very rapidly during the star's first few 
10 4 yr (the expected duration of the Class stage, see 
1993| ). It is possible that this evolution is 



driven by a rapid decrease in the infall/accretion rate on 
the central object (Bontemps et al. |1996| , Tomisaka 1996 



evolve with time from having "bullets" to not having 
them, and from being chemically active to not being so 



(see Bontemps et al. 1996 for evidence that outflow mo- 
mentum flux systematically decreases with outflow age). 
If this evolution scenario is correct, IRAS 18148 in L483 



Henriksen et al. 1997), but although Class objects do 
have stronger infall signatures than Class I sources, there 
is no clear infall trend among Class objects themselves 
(Mardones et al. 1997 ). Further study of transition objects 
like L483 is needed to understand these earliest changes of 
stellar life, and as this work has shown, the combination 
of the chemistry and kinematics of the outflow may hold 
the key to that understanding. 

5. Summary 

We have observed the L483 core and outflow in different 
mm molecular transitions and made full maps in CO(2-l), 
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CH 3 OH(2 fc -l fc ), and H 2 CO(2 2 i-ln). With these data, we 
have studied the outflow, the core, and their relation with 
the IR cometary nebula around IRAS 18148. The main 
conclusions of our work are as follows: 

1. The 12 CO emission at ambient velocities is ex- 
tremely thick with the brightness temperature expected 
for gas at 9 K, the temperature previously estimated for 
the core gas. Outside the ambient regime the 12 CO lines 
present bright wings indicating outflow material warmer 
than the ambient gas by at least a factor of 2. A simple 
model for the C 17 O(l-0) emission towards the core center 
shows that this line is optically thin and non Gaussian 
due to the presence of two velocity components. From the 
integrated C 17 O(l-0) emission we estimate a central Ho 
column density of 3 x 10 22 cm" 3 in the inner 20". 

2. The CO outflow emission is compact and slow, with 
a total length of 0.15 pc and a kinematical age of 10 4 yr. 
Lower limits to the outflow mass, momentum, and energy 
are 0.01 M Q , 0.03 M Q km s -1 , and 2 x 10 42 erg, respec- 
tively. The CO outflow is asymmetric, with a blue lobe 
having a bright spot coinciding with the NIR nebula and 
the red lobe being weaker near IRAS 18148 and having 
a relative maximum 45" from the source. This asymme- 
try suggests that the reflection nebula around the IRAS 
source may also be asymmetric and have a more promi- 
nent blue side. 

3. The CH 3 OH emission traces a dense core with no 
appreciable outflow wing contribution, although there is 
a shift in the line velocity along the direction of the out- 
flow. The CH3OH emission is centrally concentrated on 
the IRAS position, and in the central 200" (0.1 pc) it de- 
creases radially in a manner intermediate between what 
would be expected for optically thin emission with den- 
sity power laws of r" 1 and r~ 15 . The estimated mass in 
this region is 5-10 M©. No evidence for CH3OH or SiO 
abundance enhancement is found in the outflow. 

4. The H 2 CO(2 12 -l 11 ) emission is self absorbed at am- 
bient velocities, and presents spectra with brighter blue 
peak, characteristic of inward motions, toward the central 
40". With a simple model, we estimate an average infall 
speed of 0.02 km s -1 and an infall rate of 2 x 10~ 6 M 
yr" 1 . At high velocities, the H 2 CO(2i2-ln) line presents 
bright wings in the same sense as the CO outflow wings, 
indicative of outflow acceleration. Comparing the H2CO 
and CO wing intensities we find that the H2CO abun- 
dance in the outflow regime is enhanced with respect to 
the ambient regime by a factor of 20. 

5. The combination of CO, CH 3 OH, and H 2 CO data 
shows that the L483 outflow is less extreme than other 
outflows from Class objects, like L1448-C and LI 157, 
although it has some of their characteristics, such as gas 
heating and some abundance enhancement. We therefore 
suggest that the Class source at the center of the L483 
outflow is more evolved than other Class sources, and it 
is in its transition to become a Class I object. 
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Appendix A: Separation between outflow and 
ambient cloud emission 

Estimating the outflow energetics for L483 is complicated 
because its low velocity and the presence of an extended 
component. In this appendix we discuss the method we 
have applied to correct for contamination by background 
gas seen in emission. Unfortunately, the part of the out- 
flow emission absorbed by foreground gas cannot be recov- 
ered by any simple means. 

The reason it is possible to correct for the extended 
emission in L483 is because this emission seems constant 



over the flow. This can be seen in Figure A.l, where we 
present a map of spectra with the spectrum from (—40", 
40") (a representative ambient position) superposed in 
lighter shade over each map position. As the figure shows, 
positions without outflow wings have spectra with the 
same shape, suggesting that the ambient cloud superposed 
to the outflow emission contributes everywhere with a sim- 
ilar spectrum. If the outflow emission is optically thin (as 
suggested by the lack of 13 CO (1-0) wings), this extended 
component can be subtracted out, leaving outflow-only 
emission as a residual. The origin of the extended compo- 
nent is in part ambient cloud background to the outflow 
(the foreground part causes the absorption and cannot 
be corrected for) and in part background emission from 
unrelated gas from the Aquila Rift, as discussed in sec- 
tion 3.1.1. The fact that its spectrum is rather flat topped 
suggests that this emission is partly thick with a kinetic 
temperature around 10 K, like that of the ambient L483 
emission (but less extreme given its weaker 13 CO(1-0) 
emission). Given these characteristics, the gas that does 
not appear in absorption is most likely background to the 
outflow and therefore susceptible to correction. 

To avoid adding noise in the process, we have used as 
background emission the average of all non outflow posi- 
tions, and we have subtracted this spectrum to each ob- 
served position. The result seems to contain outflow emis- 
sion only, as illustrated by the total integrated emission 



map of Figure A. 2, which is very similar to the outflow 
map in Fig. || (where no background subtraction was ap- 
plied). To these background-subtracted spectra we have 
applied the standard energetics analysis (cf. Margulis & 



Lada 1985), and have ignored any contribution in the 
range Vlsr = 4-7 km s _1 , as these velocities are contami- 
nated by self absorption. For being forced to ignore these 
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Fig. A.l. Map of 12 CO(2-l) spectra towards L483 outflow (dark shade) with the spectrum from (—40", 40") super- 
posed for comparison (light shade) . Note how positions with no outflow emission show very similar spectra. This allows 
subtracting the emission from gas unrelated to outflow. For each spectrum, the velocity scale goes from Vlsr to 10 
km s _1 and the T m b scale goes from -1.5 to 10 K. 




-30 - 



Fig. A. 2. Integrated 12 CO(2-l) intensity map after back- 
ground emission subtraction (from -2 to 12 km s^ 1 ). Note 
how all residual emission appears to be from outflow (com- 
pare with Figure ^). Lower contour and step are 6 K km 
s _1 and offsets are as in Fig. [|a. 



very low velocities, our estimates are necessarily lower lim- 
its to the real outflow parameters. 
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